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Recent advances in soft-collinear effective theory with Glauber gluons have led to the development 
of a new method that gives a unified description of inclusive hadron production in reactions with 
nucleons and heavy nuclei. We show how this approach, based on the generalization of the DGLAP 
evolution equations to include final-state medium-induced parton shower corrections for large 
processes, can be combined with initial-state effects for applications to jet quenching phenomenol¬ 
ogy. We demonstrate that the traditional parton energy loss calculations can be regarded as a 
special soft-gluon emission limit of the general QCD evolution framework. We present phenomeno¬ 
logical comparison of the SCETo-based results on the suppression of inclusive charged hadron and 
neutral pion production in ^/snn = 2.76 TeV lead-lead collisions at the Large Hadron Collider to 
experimental data. We also show theoretical predictions for the upcoming y/SNN — 5.1 TeV Pb-|-Pb 
run at the LHC. 


I. INTRODUCTION 

Understanding parton shower formation and evolution is central to perturbative Quantum Chromodynamics 
(pQCD) [1]. Parton showers connect the short-distance physics of hard, large Q^, scattering with the long-distance 
non-perturbative dynamics of hadronization, govern the formations of jets [2, 3], and control the evolution of parton 
distribution and fragmentation functions through the standard Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) 
evolution equations [4-6]. High energy scattering processes naturally present a multi-scale problem, ideally suited 
to effective field theory treatment. Indeed, over the past decade many of the advances in understanding parton 
shower formation and the resummation of large logarithms that arise from ratios of energy and momentum scales in 
e+d-e”, e-bp and p-f-p have come from the well-established soft-collinear effective theory (SCET) [7-11], an effective 
theory of QCD for jet physics. Recently, SCET has been extended to describe jet propagation in matter, where 
the physical interactions with the medium are mediated via Glauber gluon exchange [12-14]. The resulting effective 
theory SCETq has been used to derive all 0{as) 1 —)■ 2 medium-induced splitting kernels [15] and discuss higher 
order 0{a^) corrections to the medium-modified jet substructure [16]. 

Such advances in the theory of in-medium parton shower formation have allowed us to overcome some of the inherent 
limitations of the traditional energy loss approach in nucleus-nucleus (A-l-A) collisions, see for example [17], and to 
unify our understanding of energetic particle and jet production in p-bp and A-bA [18]. In the current paper, we 
provide the details of the implementation of the in-medium QCD evolution-based framework and give an extended 
discussions of the connection between the energy loss approach and our new method for evaluating hadron production 
in the soft gluon emission limit. To address charged hadron and neutral pion production in Pb-bPb reactions at 
y/SNN = 2.76 TeV, we combine the medium-modified fragmentation functions with initial-state cold nuclear matter 
effects. We find that this theoretical model gives a good description of ALICE, ATLAS and CMS experimental 
measurements [19-22] of jet quenching, the attenuation in the production rate of energetic particles and jets in heavy 
ion reactions relative to the p-bp baseline scaled by the number of elementary nucleon-nucleon interactions introduced 
in [23]. While limited theoretical predictions for open heavy flavor at higher center-of-mass energies have been made 
available [24], results on light hadron production are generally absent from the literature. We take this opportunity 
to present theoretical predictions for the anticipated Pb-bPb run at — 5.1 TeV at the LHC. We note that the 

in-medium evolution approach has been previously applied to semi-inclusive deep inelastic scattering [25, 26]. 

The rest of our paper is organized as follows. In Sec. H, we discuss the strengths and limitations of the traditional 
energy loss approach to facilitate our comparison with the new QCD evolution technique. In Sec. HI, we provide 
details on how to derive the full medium-induced splitting functions and fix the virtual correction pieces through flavor 
and momentum sum rules. We demonstrate the connection between the energy loss approach and the QCD evolution 
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framework in the soft gluon approximation. In Sec. IV, we combine the final-state in-medium parton shower evolution 
with initial-state effects to provide detailed comparison to the experimental data on inclusive hadron production at the 
LHC. We conclude our paper in Sec. V. The full result for the medium-induced splitting kernels, some details about 
the phenomenological implementation and numerical evaluation of these kernels, and details about the implementation 
of cold nuclear matter effects are collected in the Appendix. 


II. THE ENERGY LOSS-BASED APPROACH 


High transverse momentum production of energetic particles in relativistic heavy ion collisions has played an essen¬ 
tial role in probing the properties of the quark-gluon plasma (QGP) at the Relativistic Heavy Ion Collider (RHIC) and 
the Large Hadron Collider (LHC). It was realized early on that as energetic quarks and gluons propagate through the 
QGP, they will interact with the medium and lose energy via collisional and radiative processes. The related attenua¬ 
tion of particle flux, commonly referred to as jet quenching [23, 27] has attracted tremendous attention from both the 
theoretical and experimental communities. There has been great progress in studying jet quenching phenomena based 
on the parton energy loss picture, ranging from inclusive light hadron and heavy meson production [19-21, 28-36] to 
di-hadron correlations [37-40] and 7 -I-hadron correlations [41-44]. 

An important step forward in understanding hard processes in the presence of a QCD medium was to realize 
that reconstructed jets will also be modified [45] and to develop a theory that describes and connects jet cross 
section and jet substructure related observables in heavy ion reactions [46]. Theoretical predictions compared to 
experimental measurements of inclusive jets [47-49], di-jets [48, 50-52], Z^/y-tagged jets [53-56] and heavy flavor 
jets production [57-59] have allowed characterization of the angular and momentum distributions of the soft medium- 
induced gluon bremsstrahlung. One approximation that is inherent in radiative energy loss-based predictions is that 
the fractional energy loss for the parent partons per emitted gluon is small and, thus, the parent partons do not 
change their identity. This limit is usually referred to as the soft gluon limit (or small-a: limit, where x = k'^/p^ 
is the light-cone momentum fraction carried away by the gluon). Large energy loss in this case proceeds through 
multiple gluon emission. Large energy loss can also result from collisional drag-like processes for the in-medium 
parton shower [60, 61] both in the weak and the strong coupling limits. 

Two approaches to parton energy loss, the High Twist (HT) approach [62, 63] and the Guylassy-Levai-Vitev (GLV) 
approach [64, 65] generalized to massive partons by Djordjevic [ 66 ], start with the picture of factorization of hard 
processes in QCD and treat the medium-induced splittings as process-dependent radiative corrections to the factorized 
expressions for hadronic, leptonic and jet observables. This makes it conceptually simpler to generalize the soft gluon 
emission limit to full parton shower treatment. We briefly discuss the energy loss approach on the example of the GLV 
solution for the medium-induced radiative spectrum. To the Hrst order in opacity [67], the distribution of inclusive 
gluon radiation is given by 
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where Cr is the quadratic Casimir in the fundamental or adjoint representation of SU(3) when the parent parton is a 
quark or a gluon, respectively. In Eq. (1) Ag is the gluon scattering length of 0(1 fm), x = k'^/p^ « u/E is the frac¬ 
tional energy of parent parton carried by the emitted gluon, is its momentum transverse to the jet axis while q^ is 

1 ^2^medium 

the transverse momentum exchange between the propagating parton and the QGP, and — —-is the normalized 

differential distribution of the elastic scattering. The details of evaluating such a double differential distribution in 
Eq. (1) are given in, e.g. Refs. [46, 47, 67]. From this distribution one can obtain the gluon number/intensity spectra, 
as well as the average gluon number/parton energy loss: 


dx J dxd?ki_ ’ 



(AE) = [ dxxE^^^^^. 
J dx 


( 2 ) 


There has also been recent interest in the scale dependence of the transport properties of matter probed in scattering 
and energy loss processes [ 68 ]. 

Despite the success of jet quenching phenomenology described above, it is now clear that qualitative advances 
in understanding QGD in the heavy ion environment are needed. Radiative energy loss calculations cannot be 
systematically improved because this approach is well defined only in the soft gluon emission limit when the parent 
parton does not change its identity. Thus, the concept of a higher order energy loss calculation is not well defined. 
Furthermore, for quarks and gluons to lose a sizable fraction of their energy would require the process to take place 
through multiple gluon emission. Precluding hard splitting processes can significantly affect the accuracy of the 
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theoretical predictions for di-hadron correlations, tagged jets and di-jets [18]. Detailed understanding of in-medium 
showers is also a prerequisite for incorporating collisional energy loss. Most importantly, the energy loss approach 
cannot incorporate the advances in understanding higher order calculations and resummation in the standard pQCD 
and SCET frameworks. 

Given the success of parton energy loss phenomenology, it is important to also validate this vast body of existing 
work and demonstrate the connection to the general QCD evolution framework in the presence of a medium. In 
Ref. [18] it was shown numerically in great detail that for sufficiently inclusive observables, such as high hadron 
production, the two approaches give practically identical results when multiple gluon emission is considered [69]. Any 
differences can be re-absorbed in a ;< 5% change of the coupling g between the jet and the medium. We demonstrate 
below that the energy loss approach is related to a special soft gluon limit solution to the DGLAP evolution equations 
with full splitting functions, when the medium-induced corrections are included. 


III. SPLITTING FUNCTIONS AND QCD EVOLUTION OF FRAGMENTATION FUNCTIONS 

In this section we review the splitting functions in both the vacuum and a QCD medium. These govern the DGLAP 
evolution equations of fragmentation functions which encode the effect of multiple gluon emissions. We will show that 
this QCD-evolution based formulation is connected to the energy loss approximation in the soft-gluon limit. Since 
the evolution equations do not make this approximation, they consistently take into account the full-a; effects in the 
medium, including flavor changing processes g ^ q,q. 


A. Splitting functions and DGLAP evolution in the vacuum 


First, we will review the leading order splitting functions in the vacuum. These are well-known and reproduced 
using the SCET collinear Lagrangian [15, 70]. The real emission contribution of the splitting kernels are calculated 
from the tree level diagrams of 1 —)■ 2 splitting. In the splitting a ^ be we define x to be the momentum fraction of 
the second final state parton c. For example, in the splitting q —)■ qg, the gluon carries the fraction x of the parent 
quark momentum. This is opposite, x 1 — a;, to traditional notation but ensures that a; ^ 1 corresponds to the 
soft gluon emission limit. With this convention, the q ^ qg splitting kernel reads: 


1 + {1-X)^ 1 ^ J_pvac, real/ N 

dxcPkx 27r2 ^ x k\ 2^“^ k\ ’ 


( 3 ) 


where (i) = q ^ qg and the subscript “vac” represents the splitting functions in the vacuum, and 
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is the corresponding leading-order splitting function. Similarly, 
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Note that the splitting functions are labeled as for a —>■ be splittings. In the energy loss approximation, 

the parent parton a emits away the soft parton c, and the parton b subsequently fragments into hadrons. 

To include the virtual contributions, we use the flavor and momentum sum rules [6] to fix the coefficients of the 
i5-function pieces in the splitting functions. The splitting functions combining real and virtual contributions take the 
following forms. 
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where we have written the real emission contributions as plus functions. The last two splitting functions do not have 
(5-function pieces since there does not exist a one-loop virtual diagram which would change the flavor of the parton. 
The flavor and momentum conservation sum rules are given by, 


[ Pqq‘'{x)dx = 0, 

Jo 

f [nT(^) + nT(^)] (1 - ^) dx = 0 , 

^0 


( 12 ) 

(13) 


/ [2nyP-(:r) + P;-(a:)] (1 -x)dx = 0, (14) 

Jo 

with n/ the number of active quark flavors. Because Pqg^’ and Pqg^’ are related hj x 1 — x, the second sum 
rule follows from the first one. Therefore we obtain 


A =0, 

P =1 dx' I^2Ca (x' + x'(l - x')) - l) - 2ny(l - x')P;“,,-(x')| 

IIC^ 2nfTR 
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Note that our plus function pieces are defined slightly differently than those in the literature. In particular, we combine 
some extra x-dependent terms into the plus functions and there is a difference between [/(x)/x]_|_ and /(x)/(x) + . In 
the small-x approximation, the real contributions to the splitting functions reduce to: 
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The flavor-diagonal splitting functions dominate, and the parton does not change its flavor in this limit. 
The fragmentation functions satisfy the DGLAP evolution equations. 
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where we have relabeled the splitting functions with |fej_| = Q and z = 1 — x. It is instructive to emphasize that the 
In Q-enhancement in the above DGLAP evolution equations can be traced back to the ^/k\ behavior in over 

a broad phase space in Q, see Eq. (3). We will comment how such a connection changes in the medium case. The 
evolved fragmentation functions can be used to describe the inclusive hadron transverse momentum spectrum in p-fp 
collisions. In the next subsection we will follow a similar procedure to describe fragmentation function evolution in 
A-l-A collisions. 


B. Medium-induced splitting functions and generalized DGLAP evolution 


The real contributions of the medium-induced splitting kernels have been computed using SGETq [15] to the first 
order in opacity. Quark-gluon splitting dressed with multiple soft interactions beyond the soft emission limit has also 
been considered [71, 72]. Eollowing the same definition as in Eq. (3), we have 
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where the subscript “med” emphasizes the medium-specific contribution, and i = q ^ qg, q —> gg, g ^ qq and 
q —>■ gq. Also, we define hi{x, fej_; /3) as the ratio of the medium and the vacuum real emission splitting functions: 


kr,l3) = 
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med, real 
Q^99 


{x,k±;l3) = 2 Ca 
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1 — X 


9^99 ' 
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Here the index f3 represents the medium properties, such as the medium size and temperature, which enter the 
medium-induced splitting function calculations. We have included the results for hpx, k±; j3) in Appendix A. 

We expect the real plus virtual medium-induced splitting functions to have the following form. 
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(28) 


P^4U^, fcx; /3) =P^4t’r"\x, kr, P). 

From momentum and flavor sum rules, the coefficients A and B, which in this case are functions of fcj_ 
properties /3, are equal to. 


(29) 

and medium 


A{kF-,P) = Q, (30) 

B{k±;P) = J dx'^^2 Ca(^x'+ x'{l-x')^ -1^ hg^gg{xPk±;/3) -2nf{l-x')P^gg{x',k±;/3)'^ . (31) 

The full splitting functions in the medium are equal to the sum of the vacuum and the medium-induced ones 

P!'^^\x,kF;P)=Prix) + Pr\^,kF;P). (32) 


Note that whereas the vacuum splitting functions are functions of the momentum fraction x alone, the medium-induced 

splitting functions depend on both x and fcx- Furthermore, P™®‘^(x, fcx; /3) = ^(x, k^; /3)+Pf^^'^’‘^{x, fcj_; ,d)-|- 

can be calculated order-by-order in the correlation between the hard scattering and the parton interactions in the 
medium, also known as opacity expansion. For medium-induced branching, only ^(x, fcj_;/?) is known beyond 

the soft gluon emission limit [15] and we use this result in the paper. In the small-x limit, the full splitting functions 
in the medium reduce to, 


Pl'p^ggix,k±;f3) = 0 , 


(^^h^^^{x,kF;P)^ , 
(^^h^^’^ix,kr,P)j , 
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(34) 

(35) 


Pg'^ggix,kj_;P) = 0, (36) 

where the function lv^^^{x, k±; j3) ^ is the same for quark and gluon splitting and can be found in Appendix B. 

Different evolution equations have been discussed in the literature, the best known examples being Balitsky-Fadin- 
Kuraev-Lipatov (BFKL) [73, 74] (evolution in x) and DGLAP (evolution in virtuality or |fcx|)- In the case of large 
processes, it has been shown [75] that the lowest order DGLAP evolution generically arises in parton showers. 
Evolution equations can be derived from the probabilistic interpretation of splitting processes inside such a shower. 


^ The superscript “sga” stands for “soft gluon approximation”. 
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This probabilistic interpretation remains the same in the presence of in-medium interactions, as long as all relevant 
interference terms for the final state interactions in the medium are included properly in the full splitting functions [15]. 
In this paper we discuss final-state parton showers, where the virtuality and ~ Q are the largest in the hard 
scattering and it decreases with branching. This is explicitly taken into account in the calculation of the relevant 
in-medium slitting kernels that we use in this study. It is worth noting that in the general equation (i.e. Eq. (32)) 
they enter as process-dependent corrections to the leading vacuum splittings. 

The full splitting functions in the medium, in Eq. (32), govern the evolution of fragmentation functions in the 
medium. Parton shower cascades in the medium behave similarly to the vacuum ones [16] and the momentum scaling 
A <C 1 of the momentum transfers from the QCD medium [14], which are integrated over, does not change the chain 
of virtuality decay in the final state. This is actually a theoretical requirement from SCETq when it was originally 
derived [12, 14]. Thus, within such a cascade picture, one might write down a generalized DGLAP-type evolution 
equation in the medium as follows. 
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A couple of comments are on order. First, since we follow the same cascade approximation, the modified DGLAP 
evolution equations in the medium have the same form as those in the vacuum in Eq. (19) and (20), with only 
the vacuum splitting functions replaced by the full splitting functions P/“". Second, such modified evolution 
equations still have the logarithmic InQ-enhancement. This is because the full medium-induced splitting functions 
are equal to the sum of the vacuum and the medium-induced ones, as in Eq. (32). The vacuum ones are always 
accompanying 1/fc^ as in Eq. (3) and thus have the logarithmic enhancement. At the same time, the medium-specific 
term, i.e., is finite in the strict k± —>■ 0 limit due to the Landau-Pomeranchuk-Migdal effect [15]. 

Eventually, since the full medium-induced splittings are the sum of the two, they thus still contain the logarithmic 
enhancement from the large phase space. This justifies the usage of a similar DGLAP-type evolution equation in 
the medium, which is true at least within the cascade approximation we are taking. Similar studies along this line 
but within different theoretical framework has been considered in, e.g.. Refs. [26, 62]. In the next subsection, we will 
see how these evolution equations encode the effect of parton energy loss from multiple soft gluon emissions. 


C. From QCD evolution to energy loss 


In the soft-gluon emission limit, the approximate solutions to the evolution equations for fragmentation functions 
are connected to the energy-loss of quarks and gluons. In this limit, the evolution equations in both the vacuum and 
the medium simplify and decouple. 
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The vacuum DGLAP evolution equations can be written explicitly as 

+ Dr/:(z,Q)ln(l-z)j, (40) 

and we expand the integrand around the end point z' = 1. Note that DJf^^{z, Q) falls steeply with increasing 2 ;, and 
we define 
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as a measure of the steepness of the fragmentation functions. To modified leading logarithmic accuracy, we neglect 
the Q-dependence in n{z) and use the one-loop running of the strong coupling as{Q) = l/( 6 oln-A^—) where 69 = 

''^QCD 

{lieA — 4n/Pfl)/(127r). The evolution equations can be solved in closed forms. 
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Following a similar procedure, we find that the solutions to the generalized DGLAP evolution equations in the 
medium are, 


D^ff{z,Q) =exp 


In ( {[n{z) - 1](1 - 2 ) - ln(l - z)} 
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This equation can be further written as 
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where 




(43) 


(44) 


(45) 

(46) 


Note that the vacuum evolution and medium correction effects factorize, as can be seen in Eqs. (43) and (44), and 
the medium causes an attenuation of the fragmentation functions. Also, both ^ and (Ng)^ depend on the scale 
Q through, 


and 

AA\ /AA\ 

For small to intermediate values of z, the suppression of the fragmentation functions (and hence the inclusive hadron 
suppression factor) is controlled by the fractional energy loss (^) and the steepness n{z) of fragmentation functions. 
On the other hand, for values of z « 1, the suppression factor becomes ~ exp(—(Ag)) which is the probability of not 
emitting a gluon. Both behaviors are the same as in the energy loss formalism. This shows the connection between 
the QCD evolution approach and the energy loss formalism in the small-x approximation. 


ax' aq/' 


^ W ■ 


(48) 


IV. INCLUSIVE HADRON SUPPRESSION IN PB+PB COLLISIONS AT THE LHC 


In this section we compare our calculations to the LHC Pb+Pb run I data (A+A), and make predictions for run II 
(A+A). We focus on single inclusive hadron production in Pb+Pb collisions of the form A + B —>■ h(pT, y)+X, where 
A and B are the incoming nuclei, and h is the hadron with transverse momentum pT and rapidity y, respectively. In 
the first part of this section, we compare with the experimental data in Pb+Pb run I at a center-of-mass energy per 
nucleon-nucleon pair = 2.76 TeV. The second part is devoted to the theoretical predictions for the Pb+Pb run 

II at y/SNN = 5.1 TeV. The nuclear modification of single inclusive hadrons in Pb+Pb collisions can be quantified 
through the nuclear modification factor Raa- 


Raa{pt) 


da^/dyd'^PT 

{Ncoi\)da^p/dyd‘^PT ’ 


(49) 


where (Acoii) is the average number of binary nucleon-nucleon collisions and is computed at a given centrality, 
and da^j^/dyd^PT and dapp/dyd^pT are the differential cross section of inclusive hadron production in A+A and 
p+p collisions, respectively. Unless otherwise specified, we focus on the single hadron production around the mid¬ 
rapidity region, and consider two centrality classes: (1) the most central collisions (0-10%) with average number of 
participants Apart = 350; (2) the mid-peripheral collisions (30-50%) with average number of participants Apart = HO. 



Experimental results are often presented in different centrality classes and, for comparison, we show data sets with 
-^part closest to our calculations. 

The invariant inclusive hadron production cross section in p+p collisions can be written as 
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while the corresponding cross section in A+A collisions reads 

1 d<yAA ^ V /" jpT, = Pt/z) pined/ 

(A^coii) dyd'^pT ^ J z"^ dyd^pT, 


(51) 


In both equations, c = {g,< 7 , g}, and D^y^^{z,Q) and D^J^^{z,Q) are the parton-to-hadron fragmentation functions 
in vacuum and QGP medium evolved between the hard and lowest scales, respectively. They are obtained from the 
vacuum and medium evolution equations, as given in Eqs. (19-20) and (37-38), respectively. For the vacuum input, 
we take KKP fragmentation functions [76], and choose the factorization, fragmentation, and renormalization scales 
Q = pt^ in the following numerical calculations. 

In Eq. (51) dac (pc = PtIz) /dyd^pr^ is the hard parton production cross section, while (pc = Pt/z) /dyd^pr^ 

is the same hard parton production cross sections calculated with isospin and initial-state cold nuclear matter (CNM) 
effects. All CNM effects considered, which include dynamical nuclear shadowing [77], Cronin effect [78, 79] and initial- 
state parton energy loss [67, 80, 81], have clear physical/dynamical origin centered around the picture of multiple 
parton scattering [82], see Appendix D for implementation. 

In p-|-A collisions there is one target nucleus to generate CNM effects through the elastic, inelastic, and coherent 
scattering of the incoming and outgoing partons. In A-l-A there are two nuclei and we expect that the cross section 
modification will be a superposition of the corresponding modifications. This does not necessarily imply that the 
cross section modification in A-l-A reactions due to CNM effects is simply twice the cross section modification in p-l-A 
reactions. The precise numerical value depends on the rapidity y of the produced particles. For example, at very 
forward rapidity the Bjorken momentum fraction of one of the hard scattering partons is a;o ~ 1, which amplifies 
the effect of cold nuclear nuclear matter energy loss but eliminates coherent power corrections. The opposite is true 
for the second parton with Bjorken momentum fraction Xf, ^ 1. It is only at midrapidity, due to the symmetry in 
the problem, that the nuclear modification due to cold nuclear matter effects in A-l-A is roughly double the nuclear 
modification expected in p-l-A collisions [30, 83]. An alternative approach would be to include CNM effects in nuclear 
parton distributions, see e.g. [84, 85]. 

In the region of interest px > 5 CeV dynamical nuclear shadowing/power corrections are negligible. Cronin 
effect can play a role up to pr — 10 CeV. The effect of initial-state cold nuclear matter (CNM) energy loss can be 
amplified near kinematic bounds [86] and is the only cold nuclear matter effect relevant to our high px phenomenology. 
While at moderate transverse momenta ALICE measurements in ^/snn = 5.02 TeV p-|-Pb collisions [87] are only 
compatible with small cold nuclear matter energy loss [88], ATLAS jet measurements [89] for the same system allow 
for considerably larger CNM effects [90]. It is, thus, important to explore at least a limited range of possibilities for 
cold nuclear matter energy loss. 


A. Comparison to LHC Pb+Pb run I data 

We are now ready to compare our calculations to the existing experimental data from the LHC Pb-|-Pb run I on 
single inclusive hadron production. We will present calculations both with and without cold nuclear matter energy loss. 
There is a trade off in the nuclear modification factor Raa of inclusive hadrons between the in-medium modification of 
fragmentation functions, driven by the coupling g between the hard partons and the QCP medium, and the strength 
of cold nuclear matter energy loss. In the region px < 50 CeV, where the experimental error bars are the smallest, 
a calculation with cold nuclear energy loss and smaller coupling g gives very similar description to the calculation 
without cold nuclear matter energy loss and slightly larger coupling g. We cannot favor one scenario versus the other 
based upon existing inclusive hadron suppression data alone. As we will show in the next sub-section, by going to 
higher px with LHC Pb-|-Pb run H the differences become more pronounced and this might provide an opportunity to 
better separate those scenarios. Of course, keeping the coupling constant g fixed, for example g = 1.9, and performing 
calculations with and without cold nuclear matter energy loss will result in slightly different magnitude of Raa, as 
we will demonstrate at the end of this subsection. 

In Fig. I, we present our results for the nuclear modification factor Raa for charged hadrons compared to data from 
Pb-|-Pb collisions at the LHC y/sjvN = 2.76 TeV measured by the ATLAS collaboration. In the top row of panels we 
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FIG. 1. The nuclear modification factor Raa for charged hadrons is calculated in central Pb+Pb collisions at the LHC 
y/SNN = 2.76 TeV, and compared with the ATLAS experimental data [21]. Top left panel: central collisions without CNM 
energy loss, and bands represent the variation of the coupling strength g = 2.0 ± 0.1. Top right panel: central collisions with 
CNM energy loss, and bands represent the variation of the coupling strength g — 1.9 ± 0.1. Middle panels are the same as the 
top panels but calculated in mid-peripheral Pb-fPb collisions. Bottom panels show the comparison of our calculations with 
ATLAS charged hadron suppression measurements in central (left), as well as mid-peripheral (right) Pb-fPb collisions. The 
dashed cyan curves are for the coupling strength g = 2.0 without CNM energy loss, while the solid magenta curves are for the 
coupling strength g = 1.9 with CNM energy loss. 


show the comparison with central (0-5%) collisions, while the middle row shows mid-peripheral (30-40%) collisions. 
For these two rows, the left panel does not include CNM energy loss, while the right panel does. The uncertainty 
bands represent the variation in the coupling strength between the jet and the medium. For the case without CNM 
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energy loss it is g = 2.0 ± 0.1, with the upper (lower) edge corresponding to g = 1.9 {g = 2.1). For the case with 
CNM energy loss it is 5 = 1.9 ± 0.1, with the upper (lower) edge corresponding to g = 1.8 {g = 2.0). 

As one can see, our calculations well reproduce the suppression of inclusive charged hadron production in 0-5% 
central Pb-|-Pb collisions at the LHC measured by ATLAS collaboration. Below px = 50 GeV both evaluations give 
comparable description of the experimental data, though the inclusion of the CNM energy loss leads to slightly larger 
curvature and, thus, better agreement with the data in the high transverse momentum region. If other components 
of the calculation have negligible uncertainties, the coupling g between the jet and the medium can be constrained 
with an accuracy of 5% and the transport properties of the medium can, thus, be extracted with 20% uncertainty as 
they scale as g'^. When comparing to the 30-40% Pb-|-Pb collisions, as measured by the ATLAS collaboration, both 
sets of results describe data equally well within the statistical and systematic uncertainties. 

To better understand the relative contribution of initial-state CNM energy loss and final-state QGP effects, we 
combine the calculations with different combinations of g and CNM energy loss in the bottom two panels of Fig. 1, for 
both central (left) and mid-peripheral (right) Pb+Pb collisions. The dashed cyan curves are the calculated Raa with 
g = 2.0 in the absence of the CNM energy loss, while the solid magenta curves are the suppression with g = 1.9 that 
includes the CNM energy loss effect. We find that at relatively low px < 25 GeV, the two curves are degenerate and 
yield practically the same suppression. However, at high px > 75 GeV, the calculations with the CNM energy loss 
effect included appear to describe the data slightly better, even though the current experimental uncertainty cannot 
definitively resolve such a difference. We further illustrate such a “degeneracy” for mid-peripheral Pb-|-Pb collisions 
in the middle panels of Fig. 1: both the calculations with g = 2.0 ± 0.1 without CNM energy loss in the left panel, 
and the one with g = 1.9 ± 0.1 but with CNM energy loss in the right panel, describe the ATLAS charged hadron 
suppression equally well. 

Similarly, we compare our calculations to ALICE charged hadron production [19] in central and mid-peripheral 
Pb-FPb collisions at the LHC ^Jsmn = 2.76 TeV in the top and middle panels of Fig. 2, respectively. The top panels 
show data sets for both 0-5% and 5-10% centrality classes. Similarly, the middle panels show both 30-40% and 40-50% 
centrality classes. These experimental centralities bracket the centralities used for our theoretical calculation and, 
therefore, we expect that theoretical calculations will fall between the respective centrality classes. This is indeed the 
case in Fig. 2. In the mid-peripheral case, the middle panels of Fig. 2, the calculations are still compatible with the 
experimental data, though the suppression measured by ALICE is slightly larger than the one seen by ATLAS. 

The relative contribution of initial-state CNM energy loss and hnal-state QGP effects is again illustrated in the 
bottom panels of Fig. 2. In the transverse momentum range covered by ALICE measurements it is not possible to 
differentiate between the curves that have g = 2.0 but do not include CNM energy loss and the curves that have 
(/ = 1.9 and include CNM energy loss. 

Last but not least, we compare our results to CMS measurements of inclusive charged hadron suppression in central 
and mid-peripheral collisions. As in the comparison with ALICE data, the top panels of Fig. 3 show both 0-5% and 
5-10% centrality classes while the middle panels of Fig. 3 show the 30-40% and 40-50% centrality classes. There is a 
one-to-one correspondence between Fig. 3 for the CMS comparison to Fig. 1, for the ATLAS comparison and Fig. 2, 
for the ALICE comparison. The comparison with CMS data supports the conclusions drawn from our comparison 
with the ATLAS and ALICE data sets. 

We now move on to compare our calculations of Raa^ for neutral pion production to the ALICE measurements in 
both central and mid-peripheral Pb-|-Pb collisions at the LHC in Fig. 4. the top panels of Fig. 4 show ALICE data 
in the 0-5% and 5-10% centrality classes, while the middle panels of Fig. 4 show data in the 40-60% centrality class. 
Surprisingly, the suppression pattern for neutral pions measured by the ALICE collaboration shows different behavior 
from the charged hadron measurements: instead of showing decreasing suppression as a function of px, the neutral 
pion Raa data from ALICE indicates a constant (or even slightly increasing) suppression over the px range. 

There is very little difference in the theoretical calculations for neutral pion production and charged hadron pro¬ 
duction as seen in all the figures and theoretical predictions for y/s^N = 5.1 TeV given in the next section. This 
is expected from isospin symmetry = (7r~^ + 7r“)/2 [76], and the fact that at high transverse momenta charged 
pion production dominates the charged hadron multiplicities. It would be quite useful if additional measurements of 
neutral pion production in Pb-|-Pb reactions can be performed to higher transverse momenta for comparison with 
inclusive charged hadron suppression. 

Finally, we note that keeping the coupling constant g fixed and performing calculations with or without cold nuclear 
matter energy loss will lead to an overall difference in the magnitude of Raa- We illustrate this in Fig. 5 for g = 1.9. In 
the left panel we show the nuclear modification factor for light hadrons and the right panel contains the corresponding 
curves for neutral pions. Solid lines include CNM energy loss and dashed lines do not. 
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FIG. 2. Same as in Fig. 1 except the comparison is made to ALICE charged hadron data for the Raa for Pb+Pb collisions at 
the LHC = 2.76 TeV. 


B. Theoretical predictions for the LHC Pb+Pb run II 

The forthcoming Pb+Pb run II at the LHC is anticipated to be performed at center-of-mass energy per nucleon- 
nucleon pair y/SNN — 5.1 TeV. The higher center-of-mass energy, combined with improvements at the LHC and 
associated detector experiments, will result in better statistics and data with higher accuracy. In this section, we 
present predictions for inclusive charged hadron and neutral pion production in both central and mid-peripheral 
Pb+Pb collisions. 

We will now outline the differences between the calculations of charged hadron and neutral pion suppression 
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FIG. 3. Same as in Fig. 1 except the comparison is made to CMS charged hadron data for the Raa for Pb+Pb collisions at 
the LHC = 2.76 TeV. 


factors at y'SNN = 2.76 TeV and ^^snn = 5.1 TeV. At higher center-of-mass energies the underlying spectra of 
energetic partons are slightly harder in the region of interest, meaning that they have a smaller negative slope. When 
convolved with the fragmentation functions in Eqs. (50) and (51), at higher smaller momentum fractions z 

in the fragmentation functions are sampled, leading to smaller suppression [18]. The physics reason for 

this effect is two-fold. On one hand, the fragmentation functions for light hadrons are falling functions of z and large 
contribution to the hadronic cross section must come from small values of z. On the other hand, the partonic cross 
sections are falling functions of and large contribution to the hadronic cross section at fixed px of the hadron must 
come from small allowed values of px^ > px- Since these are competing effects, there is a probabilistic distribution of 
z around its mean finite value. When the partonic spectra are harder they allow the z distribution to shift to smaller 
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FIG. 4. Same as in Fig. 1 except the comparison is made to ALICE neutral pions data for the Raa for Pb+Pb collisions at 
the LHC y/s^N = 2.76 TeV. 


values. Second, in-medium evolution of the fragmentation functions leads to additional softening of Dy^jJ^z^Q) and 
the effect is most pronounced at large values of z, but much less pronounced at smaller values of 2 . 

At the same time, for fixed parton or hadron transverse momentum, smaller momentum fractions Xa^ xy, are probed 
in the parton distribution functions fa/Nixa,Q), fb/Ni^b^Q), leading to an increased fraction of the hard scattered 
gluons in the final state. The strength of the leading diagonal terms in the DGLAP evolution equations Eqs. (19), 
(20) is proportional to the quadratic Casimir in the fundamental (quark) and adjoint (gluon) representations, see 
Appendix A. This will lead to larger in-medium softening of gluon fragmentation functions in comparison to quark 
fragmentation functions and, consequently, larger suppression of the cross section at higher ,ysNN when the transverse 
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FIG. 5. Comparison of calculations of the nuclear modification factor of charged hadrons (left) and neutral pions (right) with 
and without cold nuclear matter energy loss. Results for central and mid-peripheral Pb-fPb collisions are shown for g — 1.9. 


momentum is held fixed. 

Finally, the small growth of the medium density at larger center-of-mass energies, see Appendix C, will also lead 
to larger suppression. 

In terms of cold nuclear matter energy loss, we anticipate that the additional suppression effect for fixed pt hadron 
or jet production will be slightly larger at smaller center-of-mass energies [90]. The effect of cold nuclear matter energy 
loss, see Appendix D, is larger when Xa,Xb are closer to the kinematic bound. We emphasize that these competing 
effects are small and we will discuss the net result for charged hadron and neutral pion quenching below. 

We present our predictions for the LHC Pb-|-Pb run II in Fig. 6. In the top panels of Fig. 6 the predicted nuclear 
modification factor Raa for charged hadrons (left) and neutral pions (right) in central (magenta) and mid-peripheral 
(cyan) Pb-|-Pb collisions at the LHC ^/snn =5.1 TeV are plotted as a function of the hadron transverse momentum 
Pt- The bands represent a variation of the coupling strength g = 2.0±0.1 and no CNM energy loss effects. The shape 
of the nuclear modification factor is very similar to the one calculated and observed in ^/smn = 2.76 TeV reactions. 
More specifically, in the case of central Pb-|-Pb collisions and charged hadron production without CNM energy loss 
at Pt = 10 GeV the predicted charged hadron suppression is 4% larger than for LHC run 1. At pt = 100 GeV, the 
predicted suppression difference is less than 1%. 

The middle panels of Fig. 6 show a similar calculations to the top ones except that the coupling g = 1.9 ± 0.1 is 
used and CNM energy loss effects are included. Compared to the top panels, this results in a slightly slower increase 
of the Raa with pt (more quenching). The nuclear modification factor calculations once again are very similar to 
what was observed at yG)vw = 2.76 TeV. A more quantitative analysis shows that at pT = 10 GeV, the predicted 
charged hadron suppression is 3% smaller than for LHG run 1. At pt — 100 GeV, the predicted suppression is 7% 
smaller. Our results are consistent with the expectation that CNM energy loss plays a more important role at higher 
transverse momenta. 

In the bottom panels of Fig. 6 we show in more detail what role CNM energy loss may play in inclusive hadron 
production. We present scenarios: (1) the calculation with g = 2.0 in absence of CNM energy loss, (2) the calculation 
with g = 1.9 and CNM energy loss effects included, and (3) same as (2) but with 50% larger CNM energy loss. The 
motivation for discussing this last case is that inclusive jet production in p-l-A collisions can be compatible with larger 
energy loss [90]. We find that the CNM energy loss effects mainly reduce the slope of the nuclear modification factor 
Raa as a function of pt, in other words, slow down the rise in Raa as pt increases. With better precision, the future 
measurements at the LHC might be able to resolve the relative contribution from the CNM energy loss and the hot 
QCD medium effects. It should be noted, though, that the expected differences are relatively small. 


V. CONCLUSIONS 

Recent developments in SCETg, an effective theory for jet propagation in matter, have allowed for qualitative 
advances in understanding particle and jet production in heavy ion reactions and quantitative control on the system¬ 
atic uncertainties that arise from the implementation of final-state inelastic scattering processes in the QGP in jet 
quenching phenomenology. A new theoretical framework to describe inclusive particle production and suppression 
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FIG. 6. The predicted nuclear modification factor Raa for charged hadrons (left) and neutral pions (right) in central (magenta) 
and mid-peripheral (cyan) Pb-pPb collisions at the LHC ^/snn = 5.1 TeV is plotted as a function of hadron transverse 
momentum pr- In top panels the bands represent a variation of the coupling strength g = 2.0± 0.1 and we do not include CNM 
energy loss effects. In the middle panels the coupling is varied g = 1.9 ± 0.1 and we do include CNM energy loss effects. In 
the bottom we present the combined results. For the coupling strength g = 2.0, we again do not include the CNM energy loss 
effects. For the coupling strength g = 1.9, we include the usual CNM energy loss constrained at RHIC, as well as an additional 
calculation where we assume 50% larger CNM energy loss effects. 


in the heavy-ion environment was laid out in [18]. Here, we focused on the analytic connection between the new 
approach, based upon generalized DGLAP evolution equations for the fragmentation functions in dense strongly- 
interacting matter, and the concept of parton energy loss for hard processes. We found that the traditional energy 
loss phenomenology is a special soft gluon emission limit of the general QCD evolution framework. With new results 
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for the medium-evolved fragmentation functions at hand, we combined them with initial-state cold nuclear matter 
effects for applications to inclusive hadron production in A-l-A reactions and jet quenching phenomenology. 

We compared the nuclear modification factor Raa for charged hadrons measured by ATLAS, ALICE and CMS 
collaborations at y/s^N = 2.76 TeV to our predictions for central and mid-peripheral collisions. Overall, we found 
very good qualitative and quantitative agreement, which improves with the inclusion of the cold nuclear matter effects 
that make the high px rise of Raa less steep. In the case of ALICE, we also compared to the neutral pion data. 
We found that even though theoretical predictions and experimental results still mostly agree within the error bars, 
there are qualitative differences in the shape of Raa- Data appears to show suppression that is flat with pT, albeit 
in a limited transverse momentum range. At high px our theory predicts smaller suppression, essentially identical 
to that of charged hadrons. It will be very illuminating to have RaaS to much higher px in the future. From 
our comparison to data we find that the coupling between the jets in the dense QCD matter created at the LHC is 
g = 2.0 ± 0.1 in the absence of CNM energy loss and is g = 1.9 ± 0.1 when CNM energy loss effects are included. 

For the second Pb-fPb run at the LHC, expected to be at — 5.1 TeV, we presented predictions for the px 

dependence of the nuclear modification factors for charged hadrons and neutral pions in central and mid-peripheral 
collisions. We showed that even though at small and intermediate transverse momenta it is not possible to differentiate 
between phenomenological results that include CNM energy loss and a weaker coupling between the jet and the medium 
and the ones that do not include CNM energy loss but have larger coupling between the jet and the medium, such 
distinction might be possible at high px with enough statistics. 

Given the need for an accurate theory to describe reconstructed jet production in heavy-ion collisions beyond the 
energy loss approach, a logical next step will be to incorporate the full SCETg in-medium splitting functions in 
the evaluation of jet cross sections [91] and jet substructure observables, such as jet shapes [92]. We leave these 
developments for future work. 
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Appendix A: Medium-induced parton splittings 


Here, we provide the expressions for the medium-induced splittings for completeness. The calculation to Hrst order 
in opacity takes into account the contribution from the splitting induced by the interactions along the trajectory of 
the parent parton and the dominant interference with the splitting induced by the large process. The topology 
and kinematics are the same for all the splitting processes. Consequently, all results can be expressed in terms of 
universal transverse momentum vectors A±, Bx, Cx, Dx and interference phases Hi,..., H 5 , defined in [14]: 


Ax = kx, Bx = kx+ xQx, Cx = kx - {I - x)qx, Dx = kx- Qx, 


Hi — H 2 = 


Bi 


Pox{l - x) 


, Hi — H 3 — 


Pox(l -x)' 


Ho — Ho = 


Cl-Bl 

Pox(l - x) 
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p^x(l-x)’ p^x{l-x) 


(Al) 
, (A2) 


where p^ = p'^ + k~^ and the parent parton has no net transverse momentum. Recall that x = k~^ /pq and kx is the 
transverse momentum of the parton carrying momentum fraction x relative to the parent parton’s direction. 

The expressions for the continuous part of the functions read: 


dN \ 
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(A3) 


where i represents one of the four splittings i = q ^ qg, q —> gq, g -A gg, g ^ qq and ^ is the corresponding 
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vacuum splitting function. Functions hi equal: 
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Here, Xq{z), Xg{z) are the scattering lengths of quarks and gluons in the medium and {l/aei) dapp^'^''^’^/d‘^q± stands 
for normalized elastic scattering cross section of a parton in the medium. Even though this quantity varies when 
the parton is a quark or a gluon, in the high energy limit, when the t— channel dominates the elastic scattering, 
this normalized cross section does not change significantly. The symmetry oi g ^ gg, g ^ qq splitting kernels under 
a; —)■ 1 — X is most easily verified explicitly by realizing that the parton scattering cross section in the medium is 
invariant under qj_ —> —Qx- 


Appendix B: The soft gluon energy loss limit 


The energy loss phenomenology, widely used to describe jet quenching at RHIC and the LHC, is self-consistent 
only in the soft gluon limit. It is instructive to verify that in this small-x limit only two of the four medium-induced 
splitting intensities survive: 
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Note that the fractional intensity is the medium-induced splitting kernel weighed by the large lightcone momentum 
fraction of the emitted soft gluon and integrated over the available transverse momentum space. The remaining two 
splitting intensities g ^ qq, q ^ 99 are suppressed in this limit by a power of the x <C 1, Specifically, they go as 
Tb[0{x/2)\, Cf[0{x/2)]. In this limit the interference structure for all medium induced splitting functions is the same. 


Factoring the vacuum splitting function in the soft gluon approximation out according to Eq. (A3) and using the 
small X limit of the fcx unintegrated medium-induced intensities in Eq. (Bl) we obtain for both quark and gluon 
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splittings identical function /i®®®-: 
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The medium-induced splitting functions in the soft gluon approximation are given in terms of this function in Eqs. (33, 
34) above. 


Appendix C: Numerical methods and medium properties 


Even though this paper is focused on understanding the similarities and differences between the energy loss approach 
and the treatment of medium-induced parton splittings on the same footing as the vacuum branchings, we pay special 
attention to the numerical approach and the results are applicable to phenomenology. 

We implement the nuclear geometry via a standard optical Glauber model. The inelastic proton-proton scattering 
cross section can be obtained from the Particle Data Group [93] or directly from measurements. For example, the 
value we use at = 2.76 TeV is = 64 mb, well within the ALICE Icr result 62.8^4® ± 1.2 mb [94]. For 

y/'sNN = 5.1 TeV we use the value = 70 mb. In the plane transverse to the collision axis the hard scattering 
processes are distributed according to the binary collision density ~ cPNcoii/(Px±. Results relevant to the LHC 
phenomenology are calculated using full numerical evaluation of the medium-induced splitting functions. In contrast 
to jet production, the medium is distributed according to the number of participants density ~ ApartSoft 
particles that carry practically all of the energy deposited in the heavy ion collision cannot deviate a lot from such 
scaling. We take into account longitudinal Bjorken expansion since transverse expansion leads to noticeable corrections 
only in the extreme transverse velocity /3 t 1 limit [95]. In our approach all relevant finite time and finite kinematics 

integrals, such as the ones over the separation between the scattering centers Azi = Zi — Zi_i, the bremsstrahlung 
gluon phase space (fc"*", k±) and the Glauber gluon transverse momentum q_L distributions are done numerically. 

The evolving intrinsic momentum and length scales in the QGP expected to be created at the LHC are determined 
as follows: we first estimate the QGP formation times at RHIC and the LHC to be tq = 0.5 fm and tq = 0.3 fm, 
respectively. Gluons dominate the soft parton multiplicities at the LHC and their time- and position-dependent 
density can be related to charged hadron rapidity density in the Bjorken expansion model [96]: 


1 d\dNg/dy) 

T (P-X.± 
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Here, dN'^^/dy = «:Apart/2 where for central Pb-|-Pb collisions n = 8.25 at ^/s = 2.76 TeV and n = 8.6 at ^/s = 
5.1 TeV. Assuming local thermal equilibrium one finds: 


T{T,yLA_) = xj.)/16C(3) , t > tq . 


(C2) 


The Debye screening scale is given by too = (/T, recalling that we work in the approximation of a gluon-dominated 
plasma and A/ = 0. The relevant gluon mean free path is easily evaluated: \g = p with cr®® = 12)1^0?^ jrr^jg. 

In our evaluation we use an effective coupling between the jet and the medium, for example g = 2 corresponds to 
as = 0.32, but the QGP-induced bremsstrahlung is calculated with a running coupling as{Q) for the emission vertex. 
We finally note that in the medium partons acquire effective thermal mass and we include it as follows k\ —>■ -|- to|, 
in the non-asymptotic temperature regime. 


Appendix D: Implementation of cold nuclear matter effects 


To leading order in the framework of factorized perturbative QCD, in p-l-p collisions, leading parton production 
can be written as follows: 


ddc 


dyd?pT^ 


= K— ^ f ^^d?kaT fa/Nixa, k^-r) f —d^kbT fb/N{xb, klT)Hab^cd{s, t, u)5{s + i + it). (Dl) 

a,b,d 


Here, s is the center-of-mass energy squared and Hab^cd{s,t,u) are the hard-scattering coefficient functions with 
s, t, u being the usual partonic Mandelstam variables. A phenomenological A-factor can be included to account for 
higher order QGD contributions, but it cancels in the calculations of the nuclear modification factor. fa,b/N{x,k‘^) 
are the parton distribution functions with longitudinal momentum fraction x and transverse component kT- We have 
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included this /c-r-dependence in order to incorporate the Cronin effect in A+A collisions. We assume a Gaussian form 
in this variable: 




D ^utI (^t) 


fb/N{^bi kfj'p) fb/N^^a) 


o ^hT/ (^t) 


(D2) 


where fa/N{xa), fb/N{xb) are the usual collinear PDFs in a nucleon. 

In A+A collisions we take into account the following cold nuclear matter effects: 

Isospin effect. It can be easily accounted for on average in the nPDFs for a nucleus with atomic mass A and Z 
protons via: 


fa.b/ 


z 

A 





(D3) 


In Eq. (D3) fa/p{x) and fa/n{x) are the PDFs inside a proton and neutron, respectively. The PDFs in the neutron 
are related to the PDFs in the proton via isospin symmetry. 

Cronin effect. It can be modeled via multiple initial-state scatterings of the partons in cold nuclei and the corre¬ 
sponding induced parton transverse momentum broadening [14]. If the parton distribution function fa,b/Ai.Xa,bi k'^ b t) 
has a normalized Gaussian form, elastic scattering induces further IcT-broadening in the nucleus as the parton traverses 


{kl^T)AA — {kl^T)pp + ~ \ ' 


/;+ 


! OiAT. 


(D4) 


Here kb^x is the transverse momentum component for the parton prior to the hard scattering, C = ln(l [30, 97], 

and we choose 6 = 0.3 GeV“^, fff sets the soft transverse momentum transfer squared and Xg = CpICAXq are the 
related gluon and quark scattering lengths to leading order. 

Dynamical shadowing. Power-suppressed resummed coherent final-sate scattering of the struck partons leads to 
shadowing effects (suppression of the cross section in the small-a; region) [98]. The effect can be interpreted as a 
generation of dynamical parton mass in the background gluon field of the nucleus [77]. It is included via: 


Xa -t Xa 



-u J’ 


Xb Xb 





(D5) 


where Xa, Xb is the parton momentum fractions inside the target nuclei and c, d are the outgoing partons that rescatter, 
represents a characteristic scale of the multiple scattering and = CfICa^I [98]. 

Cold nuclear matter energy loss. As the partons from one nucleus undergo multiple scattering in the other nucleus 
before the large process, they can lose energy due to medium-induced gluon bremsstrahlung. This effect can be 
easily implemented as a momentum fraction shift in the PDFs 


fq,g/N{Xa) -A- fq,g/N ( , 

V-i - Eeff.q.g, 


(D6) 


Ideally, Eq. (D6) should include a convolution over the probability distribution of cold nuclear matter energy loss [86]. 
However, concurrent implementation of such distribution together with the Cronin effect and coherent power correc¬ 
tions is computationally demanding. The main effect of the fluctuations due to multiple gluon emission is an effective 
reduced fractional energy loss Ceff. ~ 0.7{AE/E). We evaluate (AE/E) as in [67] and it also depends on the medium 
parameters fff and Xg. Work is underway to better understand the physics of initial-state parton showers at lower 
fixed target energies versus coherent energy loss at much higher collider energies [99]. 

We expect that, at least as an approximation, cold nuclear matter parameters that enter the evaluation of cold 
nuclear matter effects are related. For quarks, the scale of coherent scattering per nucleon was determined to be 
^g = 0(0.1 GeV^). For the default cold nuclear matter effects used in this paper we set = 2ff^lXq = 0.08 GeV^ 
and for the 50% larger nuclear effects = 2p?/Xq = 0.12 GeV^. Inclusion of these cold nuclear matter effects gives 
(PT,) / dyd^pT,. 
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